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ABSTRACT 

We fit the Doppler profiles of the He-like triplet complexes of VII and N VI 
in the X-ray spectrum of the star £ Pup, using XMM-Newton RGS data col- 
lected over ~ 400 ks of exposure. We find that they cannot be well fit if the 
resonance and intercombination lines are constrained to have the same profile 
shape. However, a significantly better fit is achieved with a model incorporating 
the effects of resonance scattering, which causes the resonance line to become 
more symmetric than the intercombination line for a given characteristic contin- 
uum optical depth t*. We discuss the plausibility of this hypothesis, as well as 
its significance for our understanding of Doppler profiles of X-ray emission lines 
in O stars. 

Subject headings: stars: early type — star: winds, outflows — techniques: spec- 
troscopic — stars: individual (£ Pup) 



1. Introduction 



X 



High resolution X-ray spectra obtained with diffraction grating spectrometers on the 
Chandra and XMM-Newton X-ray observatories have revolutionized our understanding of 
the X-ray emission of O stars in the last five years. In the canonical pic ture, the X-rays are 
emitted in plasmas heated by shocks distributed thro ughout the wind (jCassinelli fc Swank 



1983; Corcoran et al. 1993, 1994; Hillier 



et al 



t he sh o cks are created by instabil- 



ities in the radiative driving force (e.g. iLucy &: White] Il980l ; lOwocki et al.l Il988l ; ICooper 
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1994 ; iFeldmeier et al.l 119971 ). Although some stars show anomalous X-ray emission that 



can be explained by a h ybrid mechanism involving winds channell e d by magnetic fie lds 



e.g. r Sco and 6 1 Ori C, iDonati et al.l 12002 . 12009 : ICohen et al.ll2003l : iGagne et al.l 120051 ). a 



number of "normal" O stars have X-ray spectra that are mostly consistent with the wind- 
shock paradigm (e.g. ( Pup, ( Ori, and S Ori). The papers describing the first few high 



tions dWaldron & Cassinelli 


2001 


Waldron et al. 


2004 


), but more 



2002 



these problems (IKramer et al 



2003 



Doppler profile model of lOwocki fe Cohenl (l200ll. hereafter OC01) has resolv ed many of 



Cohen et al.l 120061 ; iLeutenegger et al.l 120061 ) . The main 



outstanding problem is the relative lack of asymmetry in emission line Doppler profiles, 
which, if taken at fa ce value, would imply reductions in the literature mass- loss ra tes of an 
order of magnitude (IKramer et al.ll2003l ; ICohen et al.ll2006l ; lOwocki fc Cohenl 120061 ) . 



Although there is mounting evidence from other li nes of inquiry suggesting that the lit 



erature mass-loss rates may be systemati cally too high (IMassa et al.ll2003l ; iHillier et al.ll2003 



Bouret et al.ll2005l ; iFullerton et al.ll2006l ). there are also subtle radiative transfer effects that 
could cause emission line profiles to be more symmetric than one might naively expect . 



Two effects that have been investigated in the lite rature are porosity (IFeldmeier et al. 



2002 



Oskinova et al.ll2004j . 120061 ; lOwocki fc Cohenll2006l ) and resonance scattering (llgnace fc Gayley 



2003 



hereafter IG02). Porosity could l ower the effectiv e opac ity of the wind to X-rays, thus 
symmetrizing emission lines. However, lOskinova et al.l (120061 ) and lOwocki fc Cohenl (120061 ) 
have found that the characteristic separation scale of clumps must be very large to show 
an appreciable effect on line profile shapes, which makes it difficult to achieve a significant 
porosity effect. Resonance scattering can symmetrize Doppler profiles by favoring lateral 
over radial escape of photons; it is an intriguing possibility, but to date it has not been 
tested experimentally. 

In this paper, we present evidence for the importance of resonance scattering in some 
of the X-ray emission lines in the spectrum of the O star ( Pup. We show that the blend of 
resonance and intercombination lines of two helium-like triplets in the very high signal-to- 
noise ratio XMM-Newton Reflection Grating Spectrometer (RGS) spectrum of ( Pup cannot 
be well fit by assuming that both lines have the same profile but can be much better fit by 
assuming the profile of the resonance line is symmetrized by resonance scattering. 

This paper is organized as follows: in § [2] we discuss the reduction of over 400 ks of 
XM M-Newton RGS exposure on ( Pup; in § [3] we briefly recapitulate the results of OC01 
and ILeutenegger et al.l (120061 ) for Doppler profile modelling (§ 13.11) . and we show that the 
He-like OC01 profile model does not give a good fit to the O VII and N VI triplets of ( Pup 
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(§ 13.21) ; in §[4] we generalize the results of OC01 to include the effects of resonance scattering 
as derived in IG02 (§ 14.11) . and we fit this model to the data (§ 14.21) ; in § [5] we discuss our 
results; and in §[H]we give our conclusions. 



2. Data reduction 

The data were acquired in 11 separate pointings. The first two observations were per- 
formance verification, while the rest were calibration; they are all available in the public 
archive. The ODFs were processed with SAS versio n 7.0.0 using standard p rocedures; peri- 



ods of high background were filtered out. Only RGS (Iden Herder et al.ll200ll ) data were used 
in this paper, but EPIC data are available for most of the observations. Processing resulted 
in a co-added total of 415 ks of exposure in RGS1 and 412 ks in RGS2. The observation IDs 
used and net exposure times are given in Table [H 



RGS has random systematic wavelength scale errors with a la value of ±7 mA ( den Herder et al 



2001 ). A 7 mA shift could lead to significant systematic errors in the model parameters 
measured from a line profile. Because of this, we co-add all observations using the SAS task 
rgscombine. Assuming the systematic shifts are randomly distributed, co-adding the data 
will result in a spectrum that is almost unshifted (depending on the particular distribution 
of shifts of the individual observations), but that is broadened by 7 mA; this effect is much 
easier to account for in our analysis. We have assumed that the data do not vary intrinsi- 
cally. We have not formally verified that the data show no significant intrinsic variation, but 
upon visual inspection the data do not appear to vary more than expected from statistical 
fluctuations combined with the aforementioned random systematic errors in the wavelength 
scale. 

Spectral fitting was done with XSPE C version 12 .2.1; the line profile models are imple- 



mented as local models. The C statistic (ICashlll979l ) is used instead of x because of the 



low number of counts per bin in the wings of the profiles. 

Because of the failed CCD on RGS2, we only have RGS1 data for O VII Hea. We only 
fit RGS2 data for N VI Hea because the complex falls on a chip gap for RGS1. 

For each complex we fit, we first measured a local continuum strength from a nearby 
part of the spectrum uncontaminated by spectral features. We modeled this continuum as a 
power law with an index of 2, which is flat when plotted against wavelength. When fitting a 
line profile, we fit a combination of the local continuum (fixed to the measured value) plus 
the line profile model to the data. 
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For the N VI Hea complex, we also included emission from C VI Ly/3 at 28.4656 A, 
since the red wing of this line overlaps the blue wing of the resonance line of N VI Hea. The 
model parameters for C VI Lyj3 are assumed to be the same as for N VI, and it is assumed 
to be optically thin to resonance scattering. 

Emission lines and complexes were fit over a wavelength range of A_ < A < A+. Here 
A± = A (l ± Voo/c) ± AA, where AA is the resolution of RGS at that wavelength. A is the 
shortest wavelength in the complex for A_ and the longest wavelength for A + . 



3. Best fit He-like profile model 

3.1. The profile model 

In this section we briefly recapitulate the results of OC01 for the Doppler profile of an 
X-ray emiss ion line from an O star wi nd and the extension of these results to a He-like triplet 



complex by iLeutenegger et al.l (120061 ) . 



In the physical picture of the OC01 model, the wind is a two-component fluid; the bulk 
of the wind is relatively cool material of the order of the photospheric temperature, while a 
small fraction of the wind is at temperatures of order 1-5 MK, so that it emits X-rays. The 
cool part of the wind has some continuum opacity to X-rays and can absorb them as they 
leave the wind. 

The OC01 formalism casts the line profile in terms of a volume integral over the emis- 
sivity, attenuated by continuum absorption: 

L x = 4n [ dV Vx (fi,r)e- r ^ (1) 



where i]\(fi,r) is the emissivity at the observed wavelength A and r c (/i, r) is the continuum 
optical depth to X-rays of the wind. 

The line profile can be expressed in terms of the scaled wavelength x = (A/A — tyc/v^ = 
—Vz/voo] this gives the shift from line center in the observer's frame in units of the wind 
terminal velocity. The sign convention is such that positive x corresponds to a redshift. 

OC01 derive an expression for the line profile in terms of an integral over the inverse 
radial coordinate u = R*/r [cf. their eqn. (9)]: 



Jo w W 



(2) 

fi=—x/w(u) 
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In this equation we have used the following expressions: w(u) = v(u)/voo = (1 — u)P is the 
scaled velocity, r(x, u) is the (continuum) optical depth to X-rays emitted along a ray to the 
observer, fx(u) oc u q is the filling factor of X-ray emitting plasma, and u x = min(u , 1 — 
l^l 1 /^) is the upper limit to the integral; Uq = R*/Rq is the inverse of the minimum radius of 
X-ray emission Rq, and 1 — |x| 1///3 is a geometrical cutoff for the minimum radius that emits 
for a given value of x. The integral for L x can be evaluated numerically. 

The optical depth in this expression is derived in OC01. It is written as the product of 
the characteristic optical depth r* = kM /At^R^v^ times a dimensionless integral containing 
only terms depending on the geometry. It can be evaluated analytically for integer values 
of f3. For non-integer values of /3, the optical depth must be evaluated numerically, which 
is computationally costly, and thus not convenient in conjunction with the radial integral of 
the line profile. Because of this we assume (3 = 1 throughout this paper, which is a good 
approximation for ( Pup, as well as for O stars in general. 

The interesting free parameters of this model are the exponent of the radial dependence 
of the X-ray filling factor, q; the characteristic optical depth to X-rays of the cold plasma, 
r*; and the minimum radius for the onset of X-ray emission R . 



Leutenegger et al.l (120061 ) extend this analysis to a He-like triplet complex. The only 
difference is that the forbidden-to-intercombination line ratio has a radial dependence due 
to photoexcitation of the metastable upper level of the forbidden line: 

1 _ 1 



f 

n = - 



Kr. 



Tin 



(3) 



Here 



T + 0/0 c ~°l + 2PW{r) 
is the photoexcitation rate from the upper level of the forbidden line; it depends on the 



photospheric UV flux and scales with the geometrical dilution W(r); 0* is the photoexcitation 
rate near the photosphere, so that = 2<p 3f W(r); <p c is the critical photoexcitation rate, 
which is a parameter of the ion; and P = <f)*/4>c is a convenient dimensionless parameter 
that gives the relative strength of photoexcitation and decay to ground near the star such 
that 1Z(R*) = TZ /(1 + P) . In this paper, we u se values of P calculated from TLUSTY 



stellar atmosphere models (ILanz fc Hubeny 



Values of 7Zq are taken from lPorquet et al.l (120011 ). 



2003) as described in iLeutenegger et al.l (120061 ). 



To modify the expressions for the forbidden and intercombination line profiles to account 
for this effect, the emissivity is multiplied by the normalized line ratio: 

K{r) 



77/(7 



r}i(r) 



r)(r) 



1 + K(r) 
1 

1 + TZ(r) 



(4) 
(5) 
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3.2. Best fit model 

In this section we m odel the Doppler profiles of the O VII and N VI He-like triplets 



with the He-like profile of iLeutenegger et al 



(120061 ) described in § 13.11 The forbidden line is 



very weak for these two ions, and the intercombination line profile predicted by the model 
is not very different from the resonance line profile. The main difference in the profile of 
the resonance and intercombination lines is that the extremes of the wings are somewhat 
weaker. This is because the f/i ratio reverts to the low UV flux limit at very large radii 
(> lOOi?* for O VIII for ( Pup), so the intercombination line strength is reduced by a factor 
of a few. However, this has only a small effect on the profile shape. 

Although it is weak, the predicted strength of the forbidden line is a good check on the 
consistency of the profile model. The value of the characteristic optical depth r* can have 
a strong effect on the observed f/i ratio by setting the value of R\, the radius of optical 
depth unity. However, this effect is degenerate with the value of q, the exponent of the radial 
dependence of the X-ray filling factor. 

In Figures [1] and [21 we show the Doppler profiles of the O VII and N VI He-like triplets, 
together with the best-fit models. The best-fit parameters are given in Tables [2] and O There 
are significant residuals in both fits. The N VI triplet shows stronger residuals than O VII. 
The residuals have a systematic shape: the model predicts a greater flux than the data on 
the blue wing of the resonance line and the red wing of the intercombination line, while it 
underpredicts the data in the center of the blend. 

The systematic nature of the residuals shows that the shapes of the Doppler profiles of 
the resonance and intercombination lines are different. The residuals are consistent with the 
model resonance line being too blue and therefore too asymmetric, and the model intercom- 
bination line being too red and therefore too symmetric. 

Resonance scattering has been proposed by IG02 as an explanation for the properties 
of O star X-ray emission line Doppler profiles. If it is important, it can cause significant 
symmetrization of profiles of strong resonance lines. Because this is in qualitative agreement 
with our observations, we explore this idea further. 
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4. Best-fit model including the effects of resonance scattering 

4.1. Incorporating resonance scattering into OC01 

The wind structure that gives rise to X-ray emission is likely to be extremely com- 
plex, with a highly nonmonotonic velocity field coupled to strong density variations in all 
three dimensions. Nonetheless, even in such a medium, the overall flow pattern, to first 
approximation, should still be largely dominated by the radial acceleration of a mostly ra- 
dial, supersonic outflow. In such a medium, the radiative transfer of li ne emiss i on can be 



well modelled in terms of localized escape methods first developed by ISobolevI (119601 ) . so 
we use this as our basis for investigating the possible role of resonance scattering in X-ray 
emission lines. Our approach here generally follows that taken by IG02, with two modest 
generalizations: (1) Instead of assuming a fully optically thick line, we allow the line optical 
depth to be a free parameter that can range from the optically thin to thick limits; and (2) 
instead of assuming a constant-speed expansion, we allow for a nonzero radial acceleration. 

For a purely spherical expansion with radial velocity v and radial velocity gradient 
dv/dr, let us first define an expansion anisotropy factor, 

- rdv 1 fR\ 

a = -—-!. 6 

v or 

The Sobolev optical depth along a direction cos \x can then be written in the form 

T M = 1 2 ' ( 7 ) 

1 + (Jfl 2 

where tq is a characteristic optical depth in the lateral (/x = 0) direction [as defined further 
in eqn. (fT4"j) below]. The angle- dependent Sobolev escape probability is then given by 

pO) = 



For line photons trapped within a local Sobolev resonance layer, this represents the proba- 
bility of escape per scattering. When normalized by the angle-averaged escape probability, 

ip) = ^J_p(p)d f i, (9) 

then p(n) I (p) gives the relative angle distribution of net line emission from a localized site 
of thermal photon creation. 

To include the effects of resonance scattering on X-ray line emission, we can thus gener- 
alize the basic OC01 formalism, developed for the purely optically thin emission, by simply 
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multiplying the local emissivity [cf. eqn. ([2]) or OC01 eqn. (9)] by this normalized angle 
distribution: 

-Tc(tM,u) PW 



, , : , fx(u) 
w 6 [u) 



(p> 



(10) 

fi=—x/w(u) 



For optically thin lines (r^ <C 1), the angle emission correction recovers the simple isotropic 
scaling p(n)/ (p) = 1, while for optically thick lines (r M ^> 1), it takes the form 

Pit 1 ) = 1 + VP* /^N 

The analysis by IG02 examines this case of optically thick line emission in a constant speed 
expansion, for which dv/dr = implies a = — 1 and thus 

p(jj) 3 



<P> 2 



1-^). (12) 



To allow for a non zero radial velocity gradient, let us use here a similar /3-law velocity 
form to compute the anisotropy factor defined in eqn. flBj), 

<r=- 1. 13 

1 — M 

Note that this is independent of the wind terminal speed and for a given radius depends 
only on the velocity index /5s Q b- Taking /3s b = recovers the constant expansion case 
a — —1, while /?sob > gives a > — 1 at the lower radii (u > 0), where the wind is 
still accelerating. Overall, the parameter f3sob thus provides a convenient proxy for varying 
the relative importance of flow acceleration (compared to spherical divergence) in the local 
Sobolev escape scalings of the X-ray emission. 

Our introduction here of a separate symbol, Psob, for the velocity index relevant for 
Sobolev escape reflects the notion that in such a complex flow the local regions of X-ray 
emission need not always rigorously follow the global velocity of the bulk wind outflow. In 
particular, since the X-ray emitting material is generally too highly ionized to be directly 
driven by line-opacity, it might feasibly be better modeled as having a locally flat velocity 
gradient, dv/dr = 0, which would then be represented by a velocity index /3 Soh = instead 
of the (5 « 1 used to model the overall wind outflow. On the other hand, it might also be 
that the X-ray emitting gas is hydrodynamically so tightly coupled to the mean wind that 
even on a local scale of resonance trapping, it too exhibits a similarly positive acceleration 
that is best represented by taking = j3 = 1. In the line- fitting analysis below, we thus 
consider both these options. 



Instead of the IG02 assumption of optically thick lines, our analysis also allows general 
parameterization of the line optical thickness, as set by the overall factor r that gives 
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the Sobolev optical depth in the lateral direction /i = 0. In terms of fundamental atomic 
parameters, this is given by 

t = - — -finiCr T r. (14) 

8 r e v 

Here r e is the classical electron radius, &t is the Thomson cross section, f± is the oscillator 
strength of the transition, and Tij is the ion density. Using the steady state mass loss rate 
M = 4irpr 2 v within the assumed (3 velocity law, this can be recast in a form explicitly 
showing the dependence on inverse radius u, 

TO = T 0| * T , (15) 

w (u) 



Xr e cM ( m , , 

r 0,* = 77 <? f~ ■ ( 16 ) 



where we have defined a characteristic Sobolev optical depth 

Xr e cM ( n, 
4i^4 V P 
The factor 

Hi _ m n e _ Aj q { fx , 17 v 
p n e p /i7vm p 

gives the ratio of the ion number density to the mass density. Here Ai is the abundance of the 
element relative to hydrogen, is the ion fraction, fx is the filling factor of X-ray emitting 
plasma, and yLj^rrip is the mean mass per particle. We take this ratio to be a constant with 
radius, although in principle the ion fraction and filling factor could vary. 

In this paper, we take r * as a free parameter. For general values of the optical depth, 



the an gle-averaged escape probability cannot be evaluated analytically. Fortunately, [Castor 



(120041 . pp. 128-129) gives a computationally efficient approximation (attributed to G. Ry- 
bicki) that is accurate to ~ 1.5%, so we use this approximation to calculate (p) for finite 
values of r 0j *. 

We have implemented this X-ray emission formalism as a local model in XSPEC. The 
Sobolev optical depth has angular and radial dependence as given by eqn. (jTJ) and f|T5|) . The 
additional parameters added to the OC01 model are a switch to turn on or off completely 
optically thick scattering; the characteristic Sobolev optical depth r 0i * (used when the com- 
pletely optically thick switch is off); and the value of the velocity law exponent used in 
calculating a, ftsob- 

Figures [3] and H] compare sample results for X-ray line profiles assuming various values 
of t 0) * and /3sob- The overall trend is for higher values of r * and lower values of /3s b to give 
more symmetric profiles. The trend of lower values of /3s b to give more symmetric profiles 
is what one would expect; lowering Psob suppresses the effect of the flow acceleration in 
promoting radial photon escape, thus enhancing lateral escape an symmetrizing the profile. 
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4.2. Best-fit model including resonance scattering 

In this section we fit He-like profile models including resonance scattering to the O VII 
and N VI complexes. We fit each complex twice: once assuming /?sob = 1 and once assuming 
Psob — 0. The best-fit models are shown in Figures [6l and [HJ The best-fit parameters 
are given in Tables [2] and [3J 

The O VII profile is well fit by either value of fisob- We tested goodness of fit by 
comparing the fit statistic of 1000 Monte Carlo realizations of the model to the fit statistic 
of the data; both models are formally acceptable. The fit with /3s b = 1 is better than that 
with (3 Soh = , but only by AC = 3.8, which is about 2a for one interesting parameter. The 
fit with Psob = has a significantly smaller value of To,* than the fit with Psob = 1, as would 
be expected. The fit with /Ssob = 1 is statistically consistent with the approximation that 
the Sobolev optical depth becomes infinite. 

The N VI profile is much better fit by either model including resonance scattering than 
it is by the original model. Furthermore, the model with ftsob = gives a significantly better 
fit than the model with /3s b = 1- However, neither model is formally acceptable, and even 
the /3sob = model shows residuals of the same qualitative form as the original model, albeit 
of a much lower strength. For both models including resonance scattering, the optically thick 
approximation gives a better fit than a profile with finite Sobolev optical depth. 

To test the significance of profile broadening introduced by co-adding data with ran- 
dom systematic errors in the wavelength scale, we have also fit each best-fit model with an 
additional 7 mA Gaussian broadening. In all cases, the best-fit parameters did not change 
significantly and the fit statistics were not significantly worse. Thus we conclude that our 
analysis is not strongly affected by this broadening. 

5. Discussion 

5.1. Comparison of results 

The profile fits presented in § 13.21 clearly show that the O VII and N VI He-like triplet 
complexes in ( Pup cannot be fit by models that assume the same profile shapes for the 
resonance and intercombination lines. The profile fits presented in § 14.21 show that these 
complexes can be much better fit by a model including the effects of resonance scattering. 

However, although the O VII complex is well fit by a model including the effects of 
resonance scattering, the N VI complex shows differences in profile shape between the reso- 
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nance and intercombination line that are greater than our model can reproduce, even under 
the most generous conditions (to,* — > oo,/3s b = 0). Furthermore, one would expect the two 
complexes to show relatively similar parameters; for example, since the elemental abundance 
of nitrogen appears to be roughly twice that of oxygen, one would expect the parameter r 0j * 
to be about twice as large for the fit to N VI as it is for O VII. But a fit to the N VI profile 
with (3sob = and r 0i * ~ 10 (roughly twice the value measured for O VII) would give a 
substantially worse fit than a model with infinite Sobolev optical depth, which itself has 
significant residuals. 

The fact that the apparent discrepancy between the shapes of the resonance and inter- 
combination line profiles is much greater for N VI than for O VII implies that whatever the 
symmetrizing mechanism for the resonance line is, it is significantly stronger for N VI. There 
is no obvious explanation for this in the resonance scattering paradigm. 



5.2. Plausibility of the importance of resonance scattering 

It is worth revisiting the plausibility arguments of IG02 to confirm that one would 
expect resonance scattering to be important for these ions in the wind of £ Pup. The 
relevant quantities to estimate are the Sobolev optical depth and the ratio of the Sobolev 
length to the cooling length. 

The Sobolev length is given by (e.g. lGavlevi[l995n 



L 1 + - Vth - Vth 1 (18) 
M 1 + cr /i 2 dv j dr v/rl + a/i 2 



(19) 



The cooling length is given by 

5 kAT 
2 n e X 
as derived in IG02. 

Taking the ratio, 

L c 5 kAT v v . 9 . . . 

5 kAT 4TTfi N m p 2 w 3 (u)f x . 2 , , . 

2 A M v th u 

where we have used M = AirfiNm p n e r 2 v for a smooth wind, and added a filling factor fx 
to correct for the ratio of the density of the X-ray emitting plasma to the mean density 
expected for a smooth wind. 
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Putting in some representative numbers appropriate to ( Pup, we have 

L c , <7x w 3 (u) /VI , s 

£ M m M 6 

where M 6 is the mass-loss rate in units of 10~ 6 M Q yr -1 . We have used A = 6x 10 23 ergs s" 1 cm 3 , 
AT = 2MK, fi N = 0.6, v th = SOkms" 1 , = 2500 km s" 1 , and R* = 1.4 x 10 12 cm. 

This expression is greater than unity for lateral escape except at small radii (r < 
if the filling factor is of order unity. However, if the filling factor is significantly less than 
unity, the Sobolev approximation may not be valid. 

We now consider the expected values of the characteristic Sobolev optical depth, 

r ^_ Xr e cM f nA Xr e cM f A ^ (23) 

Putting in appropriate values, we get 

T-o,* = 120 (fi^S^l) *fxMe (24) 

We give calculations of To,*/ qifx for important lines in O star spectra in Table HI We hay e 



assumed solar abundances for all elements except C, N, and O ( jAnders fc Grevessd Il989l ). 
We assumed that the sum of CNO is equal to the solar value, with carbon being negligible 
and with nitrogen having twice the abundance of oxygen; this is an estimate based on the 
observed X-ray emission line strengths. Note that the Sobolev optical depth scales with the 
wavelength of the transition; this means that the Sobolev optical depths are significantly 
smaller for an X-ray transition than they are for a comparable UV transition. It also means 
that the longer wavelength K-shell transitions of N and O will be more strongly affected by 
resonance scattering than the shorter wavelength K-shell transitions of Ne, Mg, and Si, and 
L-shell transitions of Fe, a trend that is reinforced by the high elemental abundances of N 
and O. 

Again, if the X-ray filling factors are of order unity, the characteristic Sobolev optical 
depths for the resonance lines of N VI and O VII are large, but X-ray filling factors of 
order 10~ 3 or less are sufficient to cause the lines to become optically thin. However, the 
requirement that the Sobolev length in the lateral direction be smaller than the cooling 
length is about as stringent, so if resonance scattering is important for strong lines, the 
Sobolev approximation should also be valid. 



The high filling factors required are at odds with the simple two-component fluid picture 
of the OC01 model, since the X-ray filling factors are known to be very low. However, if we 
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take the wind to be resolved into the two components on scales of the order of the Sobolev 
length, the filling factor would just be ratio of the local density to the mean density at that 
radius. This filling factor would still likely be less than unity for the X-ray emitting plasma, 
but not as low as the X-ray filling factor for the whole wind. This conjecture is a significantly 
stronger assumption than is made in OC01. 



5.3. Impact of resonance scattering on Doppler profile model parameters 



If resonance scattering is important in Doppler profile formation in the X-ray spectra of 
O stars, it may lead to a partial reconciliation with the literature mass-loss rates. The best-fit 
models for O VII have r* = 4.1, and the best-fit model for N VI has = 3.0. If we speculate 
that somehow the resonance line of N VI is even further symmetrized than predicted by 
our model, as the residuals in our best-fit model imply, the value of r* demanded by the 
intercombination line profile residuals should be somewhat higher; a reasonable guess would 
be r* ~ 4 — 5. 



These characteristic optical depths are higher than those measured by Kramer et al. 



( 120031 ) for ( Pup by applying the model of OC01 to Doppler profiles observed with the 
Chandra HETGS; the lines studied in that paper were mostly resonance lines as well. They 
are still somewhat lower than one would expect given the literature mass-loss rates; however, 
a detailed comparison with opacity calculations and mass-loss rates remains to be done. 
Furthermore, new, sophisticated analyses of UV abso rption line profile s indicate that the 



published mass-los s rates of O star winds are too high (IMassa et al.ll2003l ; iHillier et al.ll2003 



Bouret et al.l 120051 ); the most recent systematic analysis of Galactic O stars finds that for 
at least some spectr al types, the publishe d mass-loss rates must be at least an order of 
magnitude too great (IFullerton et al.l 120061 ) . but even the more conservative revisions lower 
the mass-loss rates by a factor of a few. 

The fact that our measurements have q ~ for both O VII and N VI is an important 
additional constraint. In cases in which one would like to fit a single emission line, it is 
desirable to have as few free parameters as possible. If we can assume q = 0, we reduce the 
fitting of the profile shape to two free parameters for a line with no resonance scattering 
(Ro and r*) and three free parameters for a line with resonance scattering (the previous 
parameters in addition to To,*). Furthermore, for lines where r* is large enough to obscure 
the inner part of the wind, Rq is effectively removed as a free parameter, further reducing the 
number of free parameters to one and two for non-resonance and resonance lines, respectively. 
Thus, high signal-to-noise ratio Doppler profiles with significant continuum absorption and 
no resonance scattering may provide robust measurements of the mass-loss rates of O stars. 
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A good candidate for this is the 16.78 A line of Fe XVII, which is likely not to be optically 
thick, and which is not blended with other lines. 

5.4. Future work 

Here we give a list of issues raised by this analysis that should be addressed in future 
work. 

1. The discrepancies in the fits in this paper must be resolved. The fact that we cannot 
fit the N VI profile well is unsatisfactory. The difference between the appearance of the N VI 
complex and the O VII complex requires explanation. 

2. The effect of resonance scattering on other resonance lines in the X-ray spectrum 
should be considered. Furthermore, unless we can make concrete predictions for the impor- 
tance of resonance scattering for these lines, there may be significant fitting degeneracies 
between resonance scattering and low characteristic continuum optical depths. 

3. The effect of multiple lines on resonance scattering should be explored. Of special 
importance is the calculation of the profile of a close doublet, such as Lja. In that case, the 
splitting between the two lines is of the order of the thermal velocity of the ions. 

6. Conclusions 

We have fit Doppler profile models based on the parametrized model of OC01 to the He- 
like triplet complexes of O VII and N VI in the high signal-to-noise ratio XMM-Newton RGS 
X-ray spectrum of ( Pup. We find that the complexes cannot be well fit by models assuming 
the same shape for the resonance and intercombination lines; the predicted resonance lines 
are too blue and the predicted intercombination lines are too red. This effect is what is 
predicted qualitatively if resonance scattering is important. 

We find that models including the effects of resonance scattering give significantly better 
fits. However, there is significant disagreement between the O VII and N VI profiles in the 
degree of resonance line symmetrization that is difficult to understand in the framework of 
the resonance scattering model. Nevertheless, the general trend of the resonance scattering 
model to give more symmetrized profiles provides an interesting alternative (or supplement) 
to models that assume reduced wind attenuation due to reduced mass-loss rates and/or 
porosity. 
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Table 1. List of observations with net exposure times 





texp,RGSl 


+ b 

£exp,RGS2 


ObsID a 


(ks) 


(ks) 


0095810301 


30.6 


29.8 


0095810401 


39.7 


38.3 


0157160401 


41.5 


40.2 


0157160501 


32.8 


32.8 


0157160901 


43.4 


43.4 


0157161101 


27.0 


27.0 


0159360101 


59.2 


59.2 


0159360301 


22.0 


22.0 


0159360501 


31.5 


31.5 


0159360901 


46.6 


46.6 


0159361101 


41.1 


41.0 



^XMM-Newton Observation ID. 
b Net exposure time. 
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Table 2. Model fit parameters for O VII 



Resonance 
Scattering 


/?Sob 


Q 


T* 


u 


TO,* 


G a 


n b 


C c 


MC d 


No 
Yes 


1 



-0.21 

1 5 +0 - 06 

0.15 ±0.07 


1.6 

4 1+0-3 
+ •-■--0.4 

4.1 ±0.4 


0.62 

> 0.68 

> 0.63 


> 50 
K q+3-2 


0.91 

i i i +0.03 
± -- L - L -0.04 

-■-• uz -0.03 


6.91 
6.88 ±0.07 
6.88 ±0.07 


152.2 

85.3 

89.1 


0.578 
0.730 



a G = (/ + i)/r is assumed not to vary with radius. 

b Normalization of entire complex (r + i + f) in units of 10~ 4 photons cm -2 s _1 . 
c For 83 bins. 

d Fraction of 1000 Monte Carlo realizations of model having C less than the data. 

Note. — The first row gives the best fit for a model not including resonance scattering (i.e. the model 
of OC01 and Leutenegger et al.). The second row gives the best fit for a model including resonance 
scattering with /3g b = 1, and the last row ha s Psob = 0. We used a value of P = 1.67 x 10 4 for all 



O VII profile models (ILeutenegger et al.ll2006l ). 
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Table 3. Model fit parameters for N VI 



Resonance 




















Scattering 


Psob 


Q 


T* 


u 


TO,* 


G a 


n b 


np c 


C d 


No 




-0.34 


0.5 


0.58 




0.87 


1.562 





510.5 


Yes 


1 


-0.09 


2.1 


0.50 


thick 


1.10 


1.559 


0.87 


292.2 







0.06 


3.0 


0.48 


thick 


1.15 


1.552 


1.25 


188.4 



a G = (/ + i)/r is assumed not to vary with radius. 

b Normalization of entire N VI complex in units of 10~ 3 photons cm -2 s _1 . 
formalization of C VI Ly/3 in units of 10~ 5 photons cm -2 s _1 . 
d For 117 bins. 

Note. — The first row gives the best fit for a model not including resonance 
scattering (i.e. the model of OC01 and Leutenegger et a!.). The second row 
gives the best fit for a model including resonance scattering with /3s b = 1, and 
the last row has /3g b = 0. The C VI Ly/3 line is assumed to have the same 
values of q, r*, and Uq as the N VI triplet, and is assumed not to be affected 
by reson ance scattering. We use d a value ofP=1.01xl0 5 for all N VI profile 



models (ILeutenegger et al.l 120061 ) . 



Table 4. Expected characteristic Sobolev optical depth 





A a 




A c r 0: J qi f x d 


Ion 


Transition (A) 




(io- 3 ) 



N VI 


J" 


28 78 


6599 


9 


103 






24 QD 


1478 




90 


N VTT 

iy v 11 


uyu. 


94 78 


1 ^87 9775 e 




1 9 ^7 


VII 


r 


21.60 


0.6798 


0.45 


40 




P 


18.63 


0.1461 




7 


VIII 


Lya 


18.97 


0.1387, 0.2775 c 




7, 14 


Fe XVII 




15.01 


2.517 


0.047 


11 






15.26 


0.5970 




2.5 






16.78 


0.1064 




0.5 






17.05 


0.1229 




0.6 


Ne IX 


r 


13.45 


0.7210 


0.12 


7.0 




P 


11.55 


0.1490 




1.2 


NeX 


Lya 


12.13 


0.1382, 0.2761 c 




1.2, 2.4 


Mg XI 


r 


9.17 


0.7450 


0.038 


1.6 


Mg XII 


Lya 


8.42 


0.1386, 0.2776 c 




0.27, 0.53 


Si XIII 


r 


6.65 


0.7422 


0.036 


1.1 


Si XIV 


Lya 


6.18 


0.1386, 0.2776 c 




0.19, 0.37 



^Wavelength. 



b Oscillator str engths are from CHIANTI (IDere et al.l 119971 ; 



Landi et al.l 120061 ). 



c Assumed abundance relative to hydrogen. 

d This number is calculated using eqn. f)24p assuming a mass-loss rate 
of 10~ 6 M yr" 1 . 

e Lya is a doublet. 
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Fig. 1. — O VII triplet with best-fit OC01 He-like triplet model (not including the effects of 
resonance scattering). The top panel shows the data in black with error bars and the model 
in red. The flat red line shows the assumed continuum strength. The bottom panel shows 
the fit residuals. 
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Fig. 2.— Same as Fig. U\ but for the N VI triplet. The C VI Lya line at 28.4656 A is also 
included in the fit, as well as the other fits to the N VI triplet. 
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Fig. 3. — Comparison of the influence of different values of (3 Soh on Doppler profile shape. 
All models have q = 0, u = |, and = 5. The most asymmetric model is optically thin. 
Both of the other models use the approximation that tq^ is infinite; the more asymmetric of 
the two has fisob = 1) while the least asymmetric has Psob = 0. 
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Fig. 4. — Comparison of the influence of various values of the characteristic Sobolev optical 
depth r 0i * on Doppler profile shape. All models have q — 0, uo — |, t* = 5, and /9s b = 0. In 
order from most asymmetric to least, the models have r 0i * = 0, 1, 10, and oo. 
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Fig. 5. — O VII triplet with best-fit model assuming resonance scattering with ftsob = 1- 
Scheme is as in Fig. [TJ 
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Fig. 6. — O VII triplet with best-fit model assuming resonance scattering with ftsob = 0. 
Scheme is as in Fig. [TJ 
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Fig. 7. — N VI triplet with best-fit model assuming resonance scattering with /3s b = 1- 
Scheme is as in Fig. [TJ 
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Fig. 8. — N VI triplet with best-fit model assuming resonance scattering with /3s b = 0. 
Scheme is as in Fig. [TJ 



